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accelerating the highest energy cosmic rays, resulting from the tidal disruption 
of a star or from a disk instability. The rate and power of these flares read- 
ily explains the observed flux and density statistics of UHECRs. The photon 
bursts produced by the predicted AGN flares are discussed; they may soon be 
detectable. Observations are shown to exclude that continuous jets of powerful 
Active Galactic Nuclei are the sole source of ultrahigh energy cosmic rays; the 
stringent requirements for Gamma Ray Bursts to be the source are delineated. 

Subject headings: cosmic rays, AGN, x-rays, GRB 
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qq ■ A major challenge in theoretical High Energy Astrophysics has been to explain the 

mechanism by which ultrahigh energy cosmic rays are accelerated. The difficulty stems from 
the almost contradictory requirements of the minimum magnetic field needed to confine and 
^ ■ accelerate the cosmic rays, and the maximum magnetic field consistent w ith avoiding too 
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o5 ' much synchrotron radiation and photo-pion energy loss. It was argued by IWaxmanl (Il995l ) 



that Gamma Ray Bursts could overcome this near-contradiction - if the relativistic gamma 
factors are sufficiently high - but as we show below, UHECR acceleration by GRBs alone is 
difficult to reconcile with observational data. Sufficiently powerful AGN jets and giant lobes 
of radio galaxies could also satisfy the requirements, but too few of requisite luminosity are 
observed within the GZK distance < 100 Mpc. 

We argue here that due to tidal disruption of passing stars or episodic instabilities in 
the accretion disk, even weak AGN naturally have major flares with the required frequency 
and luminosity to account for the flux of UHE cosmic rays. First, we review the basic 
requirements for UHECR acceleration. Next, we determine the constraints from UHECR 
observations; these pose a serious challenge to GRBs, and exclude steady state AGN jets, 
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as the sole accelerators of UHECRs. Then we develop the giant AGN flare scenario and 
calculate the frequency and luminosity expected for such flares; the implied UHECR flux is 
consistent with observation. We discuss the properties of the photon bursts and show that 
GLAST and future optical surveys can test our predictions. 



2. Parameters of UHECR accelerators 

The most obvious candidates for UHECR accelerators (e.g., GRB, powerful AGN) all 
have relativistic jets. Without specializing to any particular system, we consider an out- 
flow (jet) from some central engine, with mean Lorentz factor F and with variations of the 
Lorentz factor ~ T. Due to variability in the jet velocity, one gets a turbulent flow with (at 
least) mildly relativistic shocks - as seen in the jet frame. A similar analysis applies to the 
termination shock created when such an outflow enters the ISM. We consider here only the 
acceleration of protons and we do not consider non-relativistic acceleration scenarios. 

An estimate of the parameters required to accelerate UHECRs can be obtained as 
followfl Let the characteristic size of the acceleration region in the jet frame be ~ R, with a 
characteristic magnetic field strength ~ B (also measured in the jet frame). In the jet's frame, 
we require that the shock/ turbulence accelerates a cosmic ray to energy ~ E 2 q 10 20 r _1 eV. 
To do so, the turbulent cloud must confine the cosmic ray, so 

RB>3x 10 17 T~ 1 E 20 . (1) 

(Here and below we use cgs units unless specified otherwise.) Furthermore, synchrotron 
losses during a given acceleration cycle must not exceed the energy gained in the cycle. The 
minimum energy loss in a typical acceleration cycle is ~ tcRl/c times the synchroton power 
emitted by the CR in field B, resulting in the condition 

5<r 2 £ 2 - 2 . (2) 

UHECR acceleration in turbulent shocks is accompanied by a flux of photons, which 
can be an important signature of the process. The (isotropic equivalent) total power in the 
required magnetic field (Poynting luminosity) is of order 

L ~ -cT 4 B 2 R 2 > 10 45 r 2 El erg/s. (3) 



Waxman 



(1995) 'or a clear discussion. 
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If the energy in the magnetic field, protons and electrons is in equipartition, and the energy in 
electrons is emitted through synchrotron cooling in the time it takes the shock to pass through 
the magnetic cloud, a comparable luminosity is also emitted by electrons with jet-frame 
Lorentz factors 10 3 < 7 e < 10 8 B~^; the lower bound is from assuming equipartiton with the 
mildly relativistic ions, and the upper bound is from synchrotron cooling of the electrons. 
The spectra of protons and electrons (in equipartition) are expected to be roughly flat - 
equal energy per logarithmic energy interval, leading to an approximate photon spectrunO 

uL„ > 10 44 r 2 £ 2 2 , 0.01 TB eV<hu<10T MeV, (4) 

beamed into a cone of opening angle > 

Photons with energy above the threshold e 7 « 100 MeV/ (10 11 r~ 1 E 2 o) ~ lO^-E^T erg, 
have a cross section a n ~ 10 -28 cm 2 for a photo-pion production interaction, each of which 
typically decreases the proton energy by ~10%. Therefore, the condition that photo-pion 
losses do not limit further acceleration is n^a^R < 10, where n 7 is the density of photons with 
energy above e 7 and depends on the particular environment, cooling rate, etc. We obtain an 
estimate of the maximum n 7 by assuming the total jet-frame Poynting luminosity is carried 
by photons, distributed evenly over 10 logarithmic intervals. In that case n 1 rs 10~ 3 -B 2 /e 7 , 
which with ([3]) gives a sufficient condition for avoiding excessive photo-pion energy losses: 

RB 2 < lO 17 ^. (5) 

Although the precise coefficients in (TTT[2|) depend on the details of the turbulence and shocks, 
these constraints are quite general and should be robust at the level of a factor ~few. 
Conditions (jlj) and fl^D rely in addition on the assumption of approximate equipartition of 
energy between electrons and protons and thus are less general. 

The Poynting luminosity required for UHECR acceleration in a relativistic jet ([3]) need 
not translate to a lower bound on the luminosity emitted by the jet in photons since equipar- 
tition may not be valid. Nonetheless, the total power - times the beaming factor of the 
jet ~ T -2 - does provide a lower limit on the luminosity of the accretion disk which powers 
the jet. This is not a theorem; indeed, if the disk provides a s trong large scale magnetic 



field threading the black hole, the Blanford-Znayek mechanism ( iBlandford fc Znajekl 119771 ) 
can extract energy from the BH spin rather than from the accretion flow. However such a 
low-entropy scheme would not be realized in any astrophysically realistic environment. This 
can be seen explicitly using the classical Shakura-Sunyaev solution for a disk accreting at 



2 At best, equipartition is only roughly right but not accurate in detail, as evidenced by the GRB spectrum 
being far from flat. 
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approximately the Eddington rate, to obtain the maximum magnetic field strength which 
can be contained by the pressure of the accretion disk. Assuming the B-field at the BH 
attains this field strength with optimal direction, gives an upper limit on the maximum BZ 
luminosity, Lbz ; which is ~ the Shakura-Sunyaev accretion luminosity. Therefore, 10 45 -E| is 
a robust lower limit for the bolometric luminosity produced by the accretion disk powering 
a jet responsible for UHECR acceleration; beamed photon emission from the jet itself is in 
general also expected. 



3. Observational Constraints 



A successful model for UHECR acceleration must be consistent with the following ob- 
servational constraints: 

1) Energy injection rate: Estimates of the UHECR energy injection rate per l ogarit hmic 



interval range from 0.7 — 20 x 10 44 erg Mpc 3 yr x \ the former is due to IWaxmanl (119951 ) and 



the latter from the fit of iBerezinsky! (120081 ). Thus define 



ElnE 



10 44.6±0.7p 44 g erg Mpc -3 yr -l 



(6) 



2) Number of sources: In the most detailed simulation to date, iTakami et al.l (120061 ) finds a 
best fit to AGASA data on small scale anisotropy of UHECR ar rival directions for a source 
densit y n s = 77,4 xlO -4 Mpc -3 with n.4 ~ 1. This is consistent with lThe Pierre Auger Collaboration 
( 120081 ) determination that at least 61 sources contribute to the events above 57 EeV, 20 out 
of 27 of which correlate wit hin 3.2° with objects in the IVeron-Cettv fc Veronl (120061 ) AGN 



catalog closer than 75 Mpc ( The Pierre Auger Collaboration 20071 ) - Auger08.07 below, re- 



spectively. Assuming that all the uncorrelated events originate beyond 75 Mpc and that 
the ~ 5 accidentally correlated events are representative in their distance distribution gives 
the weakest lower bound: > 40 sources within 75 Mpc, corresponding n 4 > 0.3. Note that 
n s here is the number of sources visible in UHECRs which may be larger than the number 
visible via an accompanying photon jet, since magnetic deflection may un-beam UHECRs 
relative to photons. 

3) Arrival time delay of UHECRs: The arrival time delay of a UHECR with net deflection 
S6 resulting from traversing a distance D in a random weak magnetic field is r ~ \S0 2 D/c. 
Assuming that the separations between the correlated Auger UHECRs and their candidate 
sources are due to random deflections in the extragalactic magnetic field, leads to an average 
time delay (r de i ay ) = r 5 x 10 5 yr « 10 5 yr. 
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Candidate UHECR accelerators 



The observational constraints above, combined with the results of Section [2], imply that 
conventional long-lived AGN jets (quasars, blazars, BLLacs,...) cannot be the primary site 
of UHECR acceleration. There are too few AGN with the required luminosity within 75 Mpc 
in the Auger field of vi ew and only a s mall fraction of the correlated VCV galaxies satisfy 



£boi > 10 45 £f erg s" 1 flZaw et all 120081 ). The possibility that the AGN which produced the 



observed UHECRs may no longer be so powerful does not invalidate this argument: if N is 
the number of active sources within 75 Mpc required to account for the UHECRs observed 
today, then the number of sources expected at any randomly chosen time is also N. Thus 
the probability of seeing none is miniscule si nce the expected numb er is > 40. Note that 
the time delay argument originally invoked by iFarrar fc Piranl ( 120001 ) in their high-magnetic 
deflection scenario using Cen A is valid, since it was explicitly a single-source model; that 
scenario appears now to be excluded by the few-degree observed correlation between UHECR 
and many different sources, most far from Cen A (Auger07). 



As pointed out by IWaxmanl (119951 ) , the conditions outlined in Section [2] can be met by 
a GRB, if T « 300. However the observational constraints make it difficult for GRBs to 
be the so le source of super-GZ K cosmic rays. The local (isotropic equivalent) rate of long 



GRBs is (IGuetta fc Piranl 120070 r GRB = (0.05 - 0.27) Gpc -3 yr" 1 = r 10 10 



-10±0.3 



Mpc 



-3 yr -l 

In order for GRB's to give the observed energy injection rate in UHECRs the energy per 
logarithmic interval in UHECRs produced by a typical GRB should be 10 54 6 r44.6/rio erg. 
The spectrum of UHECRs produced by an individual source is expected to be and needs to 
be rather flat. It is distributed over ln(10 9 ) ~ 20 logarithmic intervals, so the observed local 
GRB rate requires each GRB to contribute a total energy in cosmic rays of order 10 56 r W fi /rin 
erg. The typical isotropic equivalent energy of photons from a GRB is ~ fewx 10 53 erg (lAmati 



20061 ). so if GRBs are responsible for all the UHECRs, they must give a factor > 100 times 
more energy to UHECRs than to photons. Having a steeper UHECR spectrum is not a 
solution, because that would only exacerbate the difficulty of producing UHECRs at 10 20 
eV. Nor is a flatter spectrum a solution because then there would be insufficient UHECRs 
at lower energies - between the ankle and the GZK regime. Furthermore, there may be 
too few local GRBs. Denoting the beaming factors for UHECRs and photons &ucr, fry, the 
effective density of observable UHECR sources is (0.05 — 0.27)6ucR' r 5 

10- 4 Mpc" 3 . Thus 

compatibility with observation is possible only if frucR ~ 1 ^ fry ~ 10~ 2 . Additionally, for 
GRBs to be the sole source of UHECRs, the correlation with nearby AGN (Auger07) must 
be a statistical fluke or an induced correlation from the large scale structure of galaxies. 
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Giant AGN flares 



We will now argue that super-Eddington flares on super-massive black holes in (weak 
sub-Eddington) AGN are good candidates to be the source of UHECRs. The mass of an 
accretion disk near a black hole of mass M (wi t hin tw ice the last stable orbit) is ~ M Q x 



(M/(lO 7 M )) n /5, following lShakura fc Sunvaevl (1 19 73( 1 and assuming an accretion rate equal 
to 1CT 3 of Eddington and alpha parameter a = 10~ 2 . This would be a weak AGN, until a 
disk instability or tidal disruption of a passing star abruptly heats the disk. With > M Q in 
this hot disk within twice the last stable orbit, the accretion would be super-Eddington and 
one would conservatively expect flares of > 0.01 M c 2 = 2 x 10 52 erg. 

Without an accretion disk surrounding the BH, the time interval over which the stellar 
debris is accreted would be expected to be IP's to 100's of times the orbit al return time, 
which is calculated to be ~ 3 x 10 6 s (jEvans fc Kochanekill989l ; IUlmerlll999l ). leading to an 
event duration of order years and a luminosity too low in comparison to ([3]) for the event to 
accelerate 10 20 eV protons. However if the BH is surrounded by an accretion disk (even a 
thin, cold one with little accretion taking place prior to the event), we expect the tidal debris 
and the accretion disk to disrupt one another in the first passage, leading to a thick disk and 
rapid accretion, within 10's of last-stable-orbit times. Because the total amount of material 
involved in the accretion event is ~ M & , residual accretion after the event is minimal, unlike 
in the case of galaxy mergers which power AGN that last 10 7 yr or longer. 

To give a concrete example, suppose that a burst producing a 10 20 eV proton is due to 
a disk instability or to the tidal disruption of a Sun-like star, in an AGN with a 3 x 10 6 M & 
black hole. The time scale of this mass injection (t, in the AGN frame) is uncertain, ranging 
from about 10 5 s for disk instabilities to ~ 3 x 10 6 s for the rate of the tidal debris return 
(jEvans fc Kochaneklll989l ; IUlmerlll999l ). Suppose that r ~ 3 and B « 3, then to satisfy both 
(JTJ) and (jSJ) we take R ~ 3 x 10 16 . The corresponding variability time scale, RT~ l /c ~ 3 x 10 5 
s, is comparable to or less than the duration of the event, as required^. With a power of 10 45 
erg/s and duration t = 10 5 to 3 x 10 6 s, the estimated energy is 10 50 for a disk instability 
event or 3 x 10 51 erg for a tidal disruption event, comfortably compatible with the total 
energy budget in both cases. 

The stellar tidal disruption mechanism is particularly attractive because it is well- 
motivated theoretically and its predictions are relatively well-constrained. The wait-time 



3 Choices of T and B implying a variability time that is longer than the duration of the event do not 
actually satisfy ([TJ) and thus would not achieve the specified £^o- Decreasing T below 3 can lead to problems 
with the variability time, for the mechanisms considered here, while increasing T generally makes it easier to 
satisfy the requirements. From Q, larger T increases the luminosity and frequency range of photon emission. 
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per galaxy for a stellar tidal disruption by its super massive black hole, in the abs ence of 
an accretion disk, is predicted to be ~ 10 4 — 10 5 vr (jMagorrian &: Tremaind Il999l ). This 
estimate has received recent confirmation by iGezari et al.l (120071 ) who have reported the UV 
and optical observation of two tidal disruption events in otherwise inactive galaxies; the 
detailed light curves and the inferred rate of tidal disruptions in inactive galaxies are consis- 
tent with theoretical predictions. The presence of an accretion disk may increase somewhat 
the rate per galaxy of st ellar tidal dis r uption s, so we denote the wait-time between giant 
AGN flares as T 4 10 4 yr. iDonley et al.l (120021 ). comparing the ROSAT All Sky Survey and 
pointed observations, report evidence for soft x-ray flares of AGN whose estimated rate 
suggests T 4 pa 0.1 6 7 , where & 7 is the photon beaming fraction. The relationship between 
the quiescent AGN luminosity and the conditions for a giant AGN flare is uncertain, so 
we take the density of ga laxies than can host giant AGN flares to be /2 10 _2 Mpc -3 , moti- 
vated by lHao et al.l (120051 ) 's determination that the density of AGN with L[H a ] > 10 5 L Q is 
0.018 Mpc -3 . Defining the rate of giant AGN flares to be Tqaf = ^6 10~ 6 Mpc~ 3 yr _1 , obser- 
vations thus give r 6 = /2/I4 ^ 20/5 7 and an apparent number density of UHECR sources 
n 4 = 0.1 r 6 r 5 6 UCR /6 7 , comfortably above the observational lower bound n 4 > 0.3. 

Taking the cosmic ray spectrum to be flat over 20 logarithmic intervals and estimat- 
ing the total energy in UHECRs in a stellar disruption-induced flare > 0.01 M Q c 2 , gives 
cIExjcr/ dlnE = e 5 il0 51 erg, with e 5 i > 1. The predicted energy injection rate in UHECRs is 
thus also comfortably compatible with observation, reproducing (jSJ) for e 51 r 6 = 0.7r 446 . If 
disk instabilities are responsible for most UHECR-producing flares, the UHECR flux cannot 
be predicted but we can infer that the rate of instabilities needed to explain ([HD is about 
30 times higher than in the tidal disruption case, since the total UHECR energy per event 
scales with the duration (to achieve the required minimum Poynting flux) . The spectrum 
of milder disk-instabilities can indicate if this is reasonable. 



6. Photon Bursts from AGN Flares 

The UHECR burst is accompanied by a photon burst of similar duration, observation 
of which will distinguish between a disk-instability and a tidal disruption origin. Unlike for 
a GRB the accretion rate is not far above Eddington, so it is reasonable to expect emission 
from the accretion disk as in a quasar, with a total luminosity pa 10 45 erg/s. Thus an 
event capable of producing 10 20 eV protons should be accompanied by a photon burst with 
a luminosity well in excess of the peak luminosity of supernovae (~ 10 43 erg/s), with an 
isotropic component from the accretion disk and a beamed component produced in the jet 
that can only be seen if the jet happens to point toward us. If equipartition obtains, the 



8 



luminosity of the beamed component would be given by ([3]) with spectrum (J3J), but since 
equipartition may not be reliable estimator, this is a less secure prediction than the bound 
on accretion luminosity. 

If AGN bursts are the sole source of UHECRs, we can estimate pgaf - the probability 
per unit volume that some galaxy is in a flaring state, i.e., the number density of Giant AGN 
Flares - as follows. Assuming the power in cosmic rays is an efficiency factor /cr times the 
accretion power in photons, Lgaf> we have 



or about ~ 3 x 1CT 5 times the density of L* galaxies. Note that if all UHECRs are produced by 
this AGN-burst mechanism, the probability of a galaxy being in the GAF state is independent 
of the flare mechanism (tidal disruption or disk instability) since the shorter duration and 
consequent lower energy injected per disk-instability event is exactly compensated by their 
greater required frequency. SDSS has viewed of order 10 8 galaxies at least twice, but has not 
compared the results of the different observations (M. Blanton, private communication). A 
search of archival SDSS data is warranted. Future large scale surveys like PANSTARS and 
LSST will see and measure the light curve of the predicted giant AGN flares, if they are as 
strong in the optical as suggested by (jlj). 

AGN flares are in principle also observable by monitoring telescopes like Swift, Galex, 
and GLAST, although with a few exceptions only blazar-like objects have been detected at 
gamma-ray energies, so the effect of the beaming-factor 6 7 > T -2 must be included when 
estimating the rate of detectable events. GLAST observes the entire sky multiple times each 
day with the GBM (8 keV-30 MeV) and LAT (20 MeV - 300 GeV). Denoting the beamed 
SED uLy Q^-p and the 1-day sensitivity (vF u ) sen , and using (jTj), the number of observable 
flares per year is 



where the duration of the event is ta days, with ta ~ 1 — 30 for tidal disruption and stellar 
tidal disruption bursts respectively. 

Thus if AGN bursts are responsible for UHECRs, we predict that GLAST will observe 
large numbers of these bursts each year. The flare duration should be long enough in either 
the tidal disruption or disk instability case, that GLAST can issue an alert to allow detailed 
observation at other wavelengths. The spectra of individual flares and the distribution of 



Pgaf = Tqaf t 



T E inE AlnE 

fcR Lqaf 



(7) 
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maximum energies among different flares will elucidate the acceleration environment by 
comparing reality with the equipartition estimate (j3J). 



7. UHECR predictions 



Unlike conventional AGN jets or GRBs, giant AGN flares have no difficulty explaining 
the UHECR observations discussed in Section [3j The cosmic ray energy injection rate is a 
prediction of the tidal-disruption scenario and explains the observed rate ([6]) for reasonable 
parameter values. Furthermore, the effective density of visible UHECR sources is estimated 
to be well above the observational lower limit and can be larger if UHECRs are not so 
beamed as the photons from the jet. 

The giant AGN flare model makes further predictions that are supported by UHECR 
observations. The correlation between UHECRs and AGN observed by Auger is of course 
a basic prediction of this scenario, but a further prediction is that the correlating AGN - 
observed today with photons - are not especially luminous, since after a flare the AGN 
quickly returns to a weakly accreting state as discussed in Section [51 Furthermore, when the 
UHECR arrival time delay is short compared to the time since the last flare, the spectrum 
should be charact eristic of a bursting source rather than of a co ntinuous source, modulo 
GZK distortions (IWaxman fc Miralda-Escudd Il996l ; iFarrarl 120071 ) . In fact, the spectrum 
of the Ursa Major cosmic ray cluster - 4-5 events in the combined AGASA-HiRes data 
( lAbbasi et al.ll2005l ; lFarrarll2005l ) which is the only example s o far of a clu ster of events from 
a single source - favors a bursting over a continuous source (iFarrarl 120071 ). These and other 
consequences of UHECRs being produced by giant AGN flares will be discussed in greater 
detail elsewhere. 



8. Summary 

Theoretical and observational constraints on sources of ultrahigh energy protons are 
derived. Conventional, powerful AGN fail on the grounds of being too rare, and GRBs 
are only viable as the sole source of UHECRs if their energy in cosmic rays is more than 
two orders of magnitude greater than in photons. Powerful, short duration bursts of AGN, 
predicted to be induced by the tidal disruption of a star or an instability in the accretion disk, 
are shown to satisfy all of the requirements for UHECR acceleration. Numerical simulations 
of stellar tidal disruption in the presence of a thin accretion disk are needed, to test our 
prediction of efficient and rapid generation of a thick disk. The photon counterparts should 



-10- 



be observable by GLAST and in archival SDSS data or the next generation of optical surveys. 
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